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£N| ' This article reviews recent developments in quantum fluid dynamics and quantum turbulence (QT) for su- 

perfluid helium and atomic Bosc-Einstein condensates. Quantum turbulence was discovered in superfluid 4 He 
in the 1950s, but the field moved in a new direction starting around the mid 1990s. Quantum turbulence is 
comprised of quantized vortices that are definite topological defects arising from the order parameter appearing 
in Bose-Einstein condensation. Hence QT is expected to yield a simpler model of turbulence than does con- 
ventional turbulence. A general introduction to this issue and a brief review of the basic concepts are followed 
by a description of vortex lattice formation in a rotating atomic Bose-Einstein condensate, typical of quantum 
fluid dynamics. Then we discuss recent developments in QT of superfluid helium such as the energy spectra 
and dissipative mechanisms at low temperatures, QT created by vibrating structures, and the visualization of 
QT. As an application of these ideas, we end with a discussion of QT in atomic Bose-Einstein condensates. 
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1 Introduction 

o ■ 
o . 

\Q . Nature is filled with fluid flow, from small scales to large scales. Our daily life on a human scale 
constantly experiences currents of air, water, and the like. The earth sustains the ocean flows 
and atmospheric circulation, leading to an abundance of life. In space, stars and galaxies flow, 
forming many large-scale structures. Most flow in nature is actually turbulent. As the velocity 
increases, flow generally changes from laminar to turbulent. Turbulence has been a mystery for a 
long time. It has been investigated not only in basic science, including physics and mathematics 
research, but also in applied sciences, such as fluid engineering and aeronautics. As Feynman 
^ ■ said, turbulence is a problem "that is common to many fields, that is very old, and that has 
not been solved." [l[ This is chiefly because turbulence is a complicated dynamical phenomenon 
with strong nonlinearity far from an equilibrium state. 

"Turbulence" reminds one of some sketches by Leonardo da Vinci. He observed the turbulent 
flow of water and drew pictures showing that turbulence has a structure comprised of vortices 
of different sizes (Fig. 1). Vortices may be a key to understanding turbulence. However, vortices 
are not well defined for a classical fluid, and the relationship between turbulence and vortices 
remains unclear. 

Independently of these studies in classical fluid dynamics, quantum fluids such as superfluid 
helium and atomic Bose-Einstein condensates (BECs) have been investigated. These systems are 
subject to quantum restrictions; the appearance of order parameters makes the rotational motion 
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Figure 1. Sketch of turbulence by da Vinci. 



exist only in the presence of quantized vortices. A quantized vortex is a stable topological defect 
with quantized circulation. Such vortices give rise to quantum turbulence (QT). Since quantized 
vortices are well defined as elements composing a turbulent flow, QT is expected to be an easier 
system to study than classical turbulence (CT) and has a simpler model of turbulence. 

This article addresses recent developments in quantum fluid dynamics and QT in superfluid 
helium and atomic BECs. In Sec. 2, basic concepts such as superfluidity, Bose-Einstein conden- 
sation, and quantized vortices are briefly reviewed. Section 3 is devoted to quantized vortices in 
a rotating BEC, addressing typical phenomena of quantum hydrodynamics. Section 4 describes 
present important topics of QT in superfluid helium. Applying these ideas, QT in atomic BECs 
is discussed in Sec. 5. Section 6 is devoted to conclusions. 



2 Superfluidity, Bose-Einstein condensation, and quantized vortices 

Quantum fluid dynamics has been studied in connection with superfluid helium for over 50 
years. Liquid 4 He enters a superfluid state below the A point (2.17 K) through Bose-Einstein 
condensation of the atoms 0|. The characteristics of superfluidity were discovered experimentally 
in the 1930s by Kapitza et al. The hydrodynamics of superfluid helium is well described by the 
two-fluid model, in which the system consists of an inviscid superfluid (of density p s ) and a 
viscous normal fluid (of density p n ) with two independent velocity fields v s and v n . The mixing 
ratio of the two fluids depends on temperature. As the system is cooled down below the A 
point, the ratio of the superfluid component increases, and the entire fluid becomes superfluid 
below about 1 K. The Bose-condensed system exhibits the macroscopic wavefunction ^(x, t) = 
(^(x, t)\e id (*- ,t > as an order parameter. The superfluid velocity field representing the potential 
flow is given by v s = (h/m)V6 with boson mass m. Since the macroscopic wavefunction must 
be single- valued for space coordinate x, the circulation T = <f v • dt for an arbitrary closed 
loop within the fluid is quantized in terms of the quantum value n = h/m. A vortex with such 
quantized circulation is called a quantized vortex. Any rotational motion of a superfluid is only 
sustained by quantized vortices. 

A quantized vortex is a topological defect characteristic of a Bose-Einstein condensate and 
is different from a vortex in a classical viscous fluid. First, its circulation is quantized and 
conserved, in contrast to a classical vortex whose circulation can have any value and is not 
conserved. Second, a quantized vortex is a vortex of inviscid superflow. Thus, it cannot decay by 
the viscous diffusion of vorticity that occurs in a classical fluid. Third, the core of a quantized 
vortex is very thin, on the order of the coherence length, which is only a few angstroms in 
superfluid 4 He and submicron in atomic BECs. Since the vortex core is thin and does not decay 
by diffusion, it is possible to identify the position of a quantized vortex in the fluid. These 
properties make a quantized vortex more stable and definite than a classical vortex. 

Early experimental studies of superfluid turbulence focused primarily on thermal counterflow, 
in which the normal fluid and superfluid flow in opposite directions. The flow is driven by an 
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injected heat current, and it is found that the superflow becomes dissipative when the relative 
velocity between the two fluids exceeds a critical value [3|. Feynman proposed that it is a su- 
perfluid turbulent state consisting of a tangle of quantized vortices [4||. Vinen later confirmed 
Feynman's proposal experimentally by showingthat the dissipation arises from mutual friction 
between vortices and the normal flow [H, 0, 0,13] ■ The quantization of circulation was also ob- 
served by Vinen using a vibrating wire technique [§] . Subsequently, many experimental studies 
have examined superfluid turbulence (ST) in thermal counterflow systems, revealing a variety of 
physical phenomena The dynamics of quantized vortices are nonlinear and nonlocal, so it 
has not been easy to quantitatively understand the experimental results on the basis of vortex 
dynamics. A breakthrough was achieved by Schwarz, who clarified the picture of ST consisting of 
tangled vortices by a numerical simulation of the quantized vortex filament model in the thermal 



counterflow [11|, [10]. However, since the thermal counterflow has no analogy in conventional fluid 
dynamics, these experimental and numerical studies are not helpful in clarifying the relation- 
ship between ST and classical turbulence (CT). Superfluid turbulence is often called quantum 
turbulence (QT), which emphasizes the fact that it is comprised of quantized vortices. 

Comparing QT and CT reveals definite differences. Turbulence in a classical viscous fluid 
appears to be comprised of vortices, as pointed out by da Vinci. However, these vortices are 
unstable, repeatedly appearing and disappearing. Moreover, the circulation is not conserved and 
is not identical for each vortex. Quantum turbulence consists of a tangle of quantized vortices 
that have the same conserved circulation. Looking back at the history of science, reductionism, 
which tries to understand the nature of complex things by reducing them to the interactions 
of their parts, has played an extremely important role. The success of solid state physics owes 
much to reductionism. In contrast, conventional fluid physics is not reducible to elements, and 
thus does not enjoy the benefits of reductionism. However, quantum turbulence is different, 
being reduced to quantized vortices. Thus reductionism is applicable to quantum turbulence. 
Consequently, QT should lead to a simpler model of turbulence than CT. 

Based on these considerations, research into quantum fluid dynamics has opened up new direc- 
tions since the mid 1990s. One new direction has occurred in the field of low temperature physics 
by studying superfluid helium. It started with the attempt to understand the relationship be- 
tween QT and CT The energy spectrum of fully developed classical turbulence is known 
to obey the Kolmogorov law in the inertial range. Recent experimental and numerical studies 
support a Kolmogorov spectrum in QT. Following these studies, QT research on superfluid he- 
lium has moved to important topics such as the dissipation process at very low temperatures, 
QT created by vibrating structures, and visualization of QT (l5l. [l6l|. Another new direction is 
the realization of Bose-Einstein condensation in trapped atomic gases in 1995, which has stimu- 



lated intense experimental and theoretical activity 17]. As proof of the existence of superfluidity, 



quantized vortices have been created and observed in atomic BECs, and numerous efforts have 



been devoted to a number of fascinating problems 18|. Atomic BECs have several advantages 
over superfluid helium. The most important is that modern optical techniques enable one to 
directly control condensates and visualize quantized vortices. Because a series of experiments on 
BECs clearly show the properties of quantum fluid dynamics [IH, H3, 21], we shall start with 
this topic. 



3 Quantized vortices in a rotating BEC 

What happens if we rotate a cylindrical vessel with a classical viscous fluid inside? Even if the 
fluid is initially at rest, it starts to rotate and eventually reaches a steady rotation with the 
same rotational speed as the vessel. In that case, one can say that the system contains a vortex 
that mimics solid-body rotationQ A rotation of arbitrary angular velocity can be sustained 



1 Thin boundary layers of the Ekman and Stewartson types can appear in a rotating classical fluid, but we do not consider 
such complicated structures here. 
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Figure 2. Profile of (a) the condensate density and (b) the phase of the 
macroscopic wavefunction when there is a quantized vortex in a trapped BEC. 
The value of the phase varies continuously from (red) to 27r (blue). 



by a single vortex. However, this does not occur in a quantum fluid. Because of quantization 
of circulation, superfluids respond to rotation, not with a single vortex, but with a lattice of 
quantized vortices. Feynman noted that in uniform rotation with angular velocity the curl 
of the superfluid velocity is the circulation per unit area, and since the curl is 2f2, a lattice of 
quantized vortices with number density no = cux\v a /K = 2S7/k ("Feynman's rule") arranges itself 
parallel to the rotation axis 0|. Such experiments were performed for superfluid 4 He: Packard et 
al. visualized vortex lattices on the rotational axis by trapping electrons along the cores (13, H3| ■ 

This idea has also been applied to atomic BECs. Several groups have observed vortex lattices 
in rotating BECs 19, 2^, 2l|, 24]. Among them, Madison et al. directly observed nonlinear pro- 
cesses such as vortex nucleation and lattice formation in a rotating 87 Rb BEC [21]. By sudden 
application of a rotation along the trapping potential, an initially axisymmetric condensate un- 
dergoes a collective quadrupole oscillation to an elliptically deformed condensate. This oscillation 
continues for a few hundred milliseconds with gradually decreasing amplitude. Then the axial 
symmetry of the condensate is recovered and vortices enter the condensate through its surface, 
eventually settling into a lattice configuration. 

This observation has been reproduced by a simulation of the Gross-Pitaevskii (GP) equation 
for the macroscopic wavefunction ^(x, t) = |^(x, t)|e* e ( x,i ) in two-dimensional [25l.l2^| and three- 
dimensional [27(] spaces. The corresponding GP equation in a frame rotating with frequency 
fi = fiz is given by 



(i - >y)h 



di 



~V 2 + y ex (x) + <?|* (x,i)| 2 - OL 2 
Zm 



*(x,t). 



(1) 



Here V ey[ is a trapping potential, and L z = —ih(xd y — yd x ) is the angular momentum along 
the rotational axis. The interparticle potential V mt is approximated by a short-range interaction 
Vi n t — 5^( x — x'), where g = Airf^a/m is a coupling constant characterized by the s-wave 
scattering length a. The term 7 indicates phenomenological dissipation 

Figure 2 shows the profile of the condensate density |^(x, t)\ 2 and the phase 6(x, t) when there 
is a quantized vortex in a trapped BEC. The density has a hole representing the vortex core. 
The phase has a branch cut between and 2ir, and the edge of the branch cut corresponds to the 
vortex core around which the phase rotates by 2tt as the superflow circulates. One can therefore 
clearly identify the vortex both in the density and the phase. 

A typical two-dimensional numerical simulation of Eq. ([1]) [2^, [2f| for the vortex lattice for- 
mation is shown in Fig. 3, where the condensate density and the phase are displayed together. 



x We used 7 = 0.03 in this simulation [25l .[2^l . 
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The trapping potential is 



-mu 2 [(l + e x )x 2 + (1 + e y )y% 



(2) 



where u = 2ir x 219 Hz, and the parameters e x = 0.03 and e y = 0.09 describe small devia- 



tions from axisymmetry corresponding to experiments 13, 21]. We first prepare an equilibrium 
condensate trapped in a stationary potential; the size of the condensate cloud is determined 
by the Thomas- Fermi radius -Rtf- When we apply a rotation with £1 = 0.7a;, the condensate 
becomes elliptic and performs a quadrupole oscillation [Fig. 3(a)]. Then, the boundary surface 
of the condensate becomes unstable and generates ripples that propagate along the surface [Fig. 
3(b)]. As stated previously, it is possible to identify quantized vortices in the phase profile also. 
As soon as the rotation starts, many vortices appear in the low-density region outside of the 
condensate [Fig. 3(a)]. Since quantized vortices are excitations, their nucleation increases the 
energy of the system. Because of the low density in the outskirts of the condensate, however, 
their nucleation contributes little to the energy and angular momentum^] Since these vortices 
outside of the condensate are not observed in the density profile, they are called "ghost vortices" . 
Their movement toward the Thomas- Fermi surface excites ripples [Fig. 3(b)]. It is not easy for 
these ghost vortices to enter the condensate, because that would increase both the energy and 
angular momentum. Only some vortices enter the condensate cloud to become "real vortices" 
wearing the usual density profile of quantized vortices [Fig. 3(d)], eventually forming a vortex 
lattice [Fig. 3 (e) and (f)]. The number of vortices forming a lattice is given by "Feynman's rule" 
no = 20/k. The numerical results agree quantitatively with these observations. 

Note the essence of the dynamics. The initial state has no vortices in the absence of rotation. 
The final state is a vortex lattice corresponding to rotational frequency O. In order to go from 
the initial to the final state, the system makes use of as many excitations as possible, such as 
vortices, quadruple oscillation, and surface waves. We refer the readers to Ref. 26 for details. 

These experimental and theoretical results demonstrate typical behavior of quantum fluid 
dynamics in atomic BECs. 




Figure 3. Dynamics of vortex lattice formation in a rotating BEC. The figure simultaneously shows both the condensate 
density and phase. The graphs are at (a) t =63 ms, (b) 305 ms, (c) 350 ms, (d) 410 ms, (e) 450 ms, and (f) 850 ms after 
the start of the rotation. The phase varies continuously from (red) to 2-7T (blue). 



1 Actually the vortex-antivortex pairs are nucleated in the low-density region. Then the vortices parallel to the rotation are 
dragged into the Thomas-Fermi surface, while the antivortices are repelled to the outskirts. 
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4 Quantum turbulence in superfluid helium 

This section reviews research on QT since the mid 1990s. Before considering QT, we begin with 
classical fluid dynamics and statistical properties of CT [2ci ]. 



4.1 Classical turbulence 

Classical viscous fluid dynamics is described by the Navier-Stokes equation, 

JU(x,t)+v(x,t)- Vv(x,t) = -ivP(x,t) + ^V 2 v(x,t), (3) 
at p 

where v(x, t) is the velocity of the fluid, -P(x, t) is the pressure, p is the density of the fluid, 
and v is the kinematic viscosity. The flow can be characterized by the ratio of the second term 
on the left-hand side of Eq. ([3]), referred to as the inertial term, and the second term on the 
right-hand side, called the viscous term. This ratio defines the Reynolds number R = vD/u, 
where v and D are a characteristic velocity and length scale, respectively. When v increases and 
the Reynolds number exceeds a critical value, the system changes from a laminar to a turbulent 
state, in which the flow is complicated and contains eddies. 



Such turbulent flow shows characteristic statistical behavior [29j, [3CJ]. We assume a steady 
state of fully developed turbulence of an incompressible classical fluid. Energy is injected into 
the fluid at a rate e in an energy-containing range. In an inertial range, this energy is transferred 
to smaller length scales without dissipation. In this range, the system is locally homogeneous 
and isotropic, leading to energy spectral statistics described by the Kolmogorov law, 

E(k) = C£ 2/3 k~ 5/3 . (4) 

The energy spectrum E(k) is defined by E = J dkE(k), where E is the kinetic energy per unit 
mass and k is the wavenumber from the Fourier transform of the velocity field. The spectrum 
of Eq. (jlj) is derived by assuming that E(k) is locally determined only by the energy flux e 
and by k. The energy transferred to smaller scales is dissipated at the Kolmogorov wavenumber 
kx = (e/V 3 ) 1 / 4 in an energy-dissipative range via the viscosity of the fluid at a dissipation rate 
e in Eq. (|4|), which is equal to the energy flux II in the inertial range. The Kolmogorov constant 
C is a dimensionless parameter of order unity. The Kolmogorov spectrum has been confirmed 
experimentally and numerically in turbulence at high Reynolds numbers. The inertial range is 
thought to be sustained by a self-similar Richardson cascade in which large eddies break up 
into smaller eddies through vortex reconnections. In CT, however, the Richardson cascade is not 
completely understood because it is impossible to definitely identify each eddy. 



4.2 Dynamics of quantized vortices 

As described in Sec. 2, most experimental studies on superfluid turbulence have examined ther- 
mal counterflow. However, the nonlinear and nonlocal dynamics of vortices delayed pro gres s in 



the microscopic understanding of the vortex tangle. Schwarz overcame these difficulties [Hi, LL2J 
by developing a direct numerical simulation of vortex dynamics connected with dynamical scaling 
analysis, enabling the calculation of physical quantities such as the vortex line density, anisotropic 
parameters, and mutual friction force. The observable quantities calculated by Schwarz agree 
with typical experimental results. 

Two formulations are generally available for studying the dynamics of quantized vortices. One 
is the vortex filament model and the other is the Gross-Pitaevskii (GP) model. 
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4-2.1 Vortex filament model 

As discussed in Sec. 1, a quantized vortex has quantized circulation. The vortex core is ex- 
tremely thin, usually much smaller than other characteristic length scales of the vortex motion. 
These properties allow a quantized vortex to be represented as a vortex filament. In classical 
fluid dynamics, the vortex filament model is an idealization. However, the vortex filament model 
is accurate and realistic for a quantized vortex in superfluid helium. 

The vortex filament model represents a quantized vortex as a filament passing through the 
fluid, having a definite direction corresponding to its vorticity. Except for the thin core region, 
the superflow velocity field has a classically well-defined meaning and can be described by ideal 



fluid dynamics [ll|, [3l|. The velocity at point r due to a filament is given by the Biot-Savart 
expression, 

k f (si — r) x dsi 



sV ; AnJc |si — r| 3 V ' 

where k is the quantum of circulation. The vortex moves with the superfluid velocity. At finite 
temperatures, mutual friction occurs between the vortex core and the normal flow. 



A numerical simulation method on this model has been described in detail elsewhere 11 



12|, |3l|]. A vortex filament is represented by a single string of points separated by distance A£. 
The vortex configuration at a given time determines the velocity field in the fluid, thus causing 
the vortex filaments to move. Vortex reconnection needs to be included when simulating vortex 
dynamics. A numerical study of a classical fluid shows that close interaction of two vortices 
leads to their reconnection, chiefly because of the viscous diffusion of vorticity. Schwarz assumed 
that two vortex filaments reconnect when they come within a critical distance of each other 
and showed that statistical quantities such as the vortex line density are not sensitive to how 



the reconnections occur Even after Schwarz's study, it remained unclear as to whether 

quantized vortices can actually reconnect. However, Koplik and Levine directly solved the GP 
equation to show that two closely quantized vortices reconnect even in an inviscid superfluid 
[32]. More recent simulations have shown that reconnections are accompanied by the emission 



of sound waves having wavelengths on the order of the healing length [33j, \34 

Schwarz numerically studied how remnant vortices develop into a statistical steady vortex 

tangle under thermal counterflow fl2]. The tangle is self-sustained by the competition between 

the excitation due to the applied flow and the dissipation by mutual friction. The numerical 

results were quantitatively consistent with experimental results, confirming Feynman's original 

picture that QT is a tangle of quantized vortices. 
The following two quantities are characteristic of a vortex tangle. The line density L is defined 

as the total length of vortex cores per unit volume. The mean spacing t between vortices is given 

by^L- 1 /?. 

4-2.2 The Gross-Pitaevskii (GP) model 

In a weakly interacting Bose system, the macroscopic wavefunction ^(x, t) appears as the 
order parameter of the Bose-Einstein condensate, obeying the Gross-Pitaevskii (GP) equation 



121, 



ih ^T 1 = ("l!l v2 + 5l * (x ' t)|2 " ") * (x ' t] - (6) 

Writing ^ = 1^1 exp(z#), the squared amplitude 1^1 2 is the condensate density and the gradient 
of the phase 6 gives the superfluid velocity v s = (h/m)V0, corresponding to frictionless flow of 
the condensate. This relation causes quantized vortices to appear with quantized circulation. The 
only characteristic scale of the GP model is the coherence length defined by £ = h/(^2mg \^\), 
which determines the vortex core size. The GP model can explain not only the vortex dynamics 
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but also phenomena related to vortex cores, such as reconnection and nucleation. However, the 
GP equation is not applicable to superfluid 4 He, which is not a weakly interacting Bose system. 
But it is applicable to Bose-Einstein condensation of a dilute atomic Bose gas 17}] , as discussed 
in Sec. 3. 



4.3 Modern research trends 

Most older experimental studies on QT were devoted to thermal counterflow. Since this flow 
has no classical analog, these studies did not significantly contribute to the understanding of the 
relation between CT and QT. Since the mid 1990s, important experimental studies have been 
published on QT in the absence of thermal counterflow, differing significantly from previous 
studies [H, [13, llB] - 



The first important contribution was made by Maurer and Tabeling [35(, who confirmed ex- 
perimentally the Kolmogorov spectrum in superfluid 4 He. A turbulent flow was produced in a 
cylinder by driving two counter-rotating disks. The authors measured the local pressure fluctua- 
tions to obtain the energy spectrum. The experiments were made at three different temperatures: 
2.3 K, 2.08 K, and 1.4 K. Both above and below the A point, the Kolmogorov spectrum was 
confirmed. 

Next was a series of experiments on g rid turbulence performed for superfluid 4 He above 1 
K by the Oregon group [H, H3, [H,j39]. Flow through a grid is usually used for generating 
turbulence in classical fluid dynamics [2a]. At sufficient distance behind the grid, the flow displays 
homogeneous isotropic turbulence. This method has also been applied to superfluid helium. In 
the Oregon experiments, the helium was contained in a channel along which a grid was pulled at 
constant velocity. A pair of second-sound transducers was implanted in the walls of the channel 
to observe vortex tangles. In combining their observations with the decay of the turbulence, the 
authors had to make several assumptions. The analysis is too complicated to be described in 
detail here. The important point is that the coupling between the superfluid and the normal 
fluid by mutual friction makes a quasiclassical flow appear at length scales much larger than 
I and causes the fluid to behave like a one-component fluid [40]. The line density is found to 
decay as i~ 3 / 2 . A simple analysis shows that the i -3 / 2 decay can be derived from a quasiclassical 
model with a Kolmogorov spectrum 



Subsequently many experimental, theoretical, and numerical works appeared. The present 
areas of study of QT are [l6| : the energy spectra and the dissipation mechanism at zero 
temperature; QT created by vibrating structures; and visualization of QT. 

4-3.1 Energy spectra and dissipation at zero temperature 

The pioneering experiments by Maurer et al. and by the Oregon group were performed at rela- 
tively high temperatures where the normal fluid component was significant, and so the situation 
is complicated because of the two-fluid behavior. Here we confine ourselves to the case of pur e 
superfluid at zero temperature. What happens to QT at zero temperature is not trivial [15l ]. 
The first problem is the nature of the energy spectrum of the turbulence for the pure superfluid 
component. The second problem is the dissipation in this system. Since there is no normal fluid 
component, any dissipative mechanism does not work at large scales. However, some dissipative 
mechanism should operate at small length scales. One possibility is acoustic emission during 
vortex reconnections. In classical fluid dynamics, vortex reconnections cause acoustic emission. 
In quantum fluids, numerical simulations of the GP model show acoustic emission at every recon- 
nection event [33l ]. However, this mechanism is thought not to be important in superfluid helium 
because of its short coherence lengthQ Another possible mechanism is the radiation of sound 
(phonons) by the oscillatory motion of the vortex cores. The third problem is how energy is 



1 Acoustic emission may be more efficient in atomic BECs in which the coherence length is generally not much shorter than 
the system size. 
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transferred from large to small scales at which the dissipative mechanism operates. If the energy 
spectrum obeys the Kolmogorov law, the energy should be transferred through a Richardson 
cascade of quantized vortices in contrast to the classical case. Vortex reconnections sustaining a 
Richardson cascade should be less effective at size scales shorter than the mean vortex spacing 

L 

No experimental studies have directly addressed the energy spectra of QT at zero temperature, 
although three numerical studies have been made. The first study was done by Nore et al. using 



the GP model [4l|, |42|] . They numerically solved the GP equation starting from Taylor-Green 
vortices and followed their time development. The quantized vortices became tangled and the 
energy spectra of the incompressible kinetic energy appeared to obey the Kolmogorov law for 
a short period, though eventually deviated from it. The second study was done by the vortex 
filament model [43)]. Araki et al. made a vortex tangle arising from Taylor-Green vortices and 
obtained an energy spectrum consistent with the Kolmogorov law. The third used the modified 
GP model by Kobayashi and Tsubota 44 , 51| . 

The Kolmogorov spectra were confirmed for both decaying [3] and steady [45| QT by the 
modified GP model. The normalized GP equation is 

i|U(x, t) = [- V 2 - M + i/|$(x, t)\ 2 Mx, t), (7) 

which determines the dynamics of the macroscopic wavefunction 3>(x, t) = /(x, t) exp[i0(x, £)]. 
The hydrodynamics in the GP model are compressible. The total number of condensate particles 
is N = f c£x|<&(x, i)| 2 and the total energy is 

E(t) = 1 J dx $*(x, t) [ - V 2 + |/(x, t) 2 ] $(x, t), (8) 

as represented by the sum of the interaction energy Ei nt (t), the quantum energy E q (t), and the 
kinetic energy E kin (t) 0,52], 

Eint(t) = 2^/ d*f(x,t)\E q (t) = ^Jd* [Vf(x,t)] 2 ,E km (t) = 1 J dx. ]/(x, t)V0(x, t)] 2 (.9) 

The kinetic energy is further divided into a compressible part Ef. in (t) due to compressible excita- 
tions and an incompressible part E l kin (t) due to vortices. The Kolmogorov spectrum is expected 
for El m (t). 

The failure to obtain a Kolmogorov law in the pure GP model [4ll . 1421 ] is attributable to 
the following. The simulations showed that E l kin (t) decreases and E? in (t) increases while the 
total energy E(t) is conserved because many compressible excitations are created through vortex 



reconnections |33l. 1341] and disturb the Richardson cascade of quantized vortices. Kobayashi et al. 
overcame these difficulties and obtained a Kolmogorov spectrum in QT that revealed an energy 
cascade 44], 4^] • They made numerical calculations of the Fourier-transformed GP equation with 



dissipation, 

7(k)) |U(k, t) = [k 2 - /i(t)]*(k, t) + ^_ $(ki, t)$*(k 2j t) x $(k - ki + k 2 , t). (10) 

ki,k 2 

Here $(k, t) is the spatial Fourier component of ^(x, t) and V is the system volume. The healing 
length is £ = l/\$\y/g. The dissipation has the form 7(k) = 7o^(fc — 27r/£) with step function 
9, which only dissipates excitations smaller than £. This form of dissipation can be justified by 
a coupled analysis of the GP equation and the Bogoliubov-de Gennes equations for thermal 
excitations 146 
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Kobayashi et al. confirmed the Kolmogorov spectra for decaying turbulence 441 ] . To obtain a 
turbulent state, they started the calculation from an initial configuration in which the density 
was uniform and the phase of the wavefunction had a random spatial distribution. The initial 
state was dynamically unstable and soon developed turbulence with many vortex loops. The 
spectrum E l kin (k,t) was then found to obey the Kolmogorov law. 

A more elaborate analysis was made for steady QT by introducing energy injection at large 



length scales in addition to energy dissipation at small scales 45)]. Energy injection at large 
scales was effected by moving a random potential V(x, t) that had a characteristic spatial scale 
Xq. Once steady QT occurs as a balance between energy injection and dissipation, it has an 
energy-containing range k < 2tt/Xq, an inertial range 2tt/Xq < k < 2tt/^, and an energy- 
dissipative range 2tt/£ < k. A typical simulation of steady turbulence was performed for Xq = 4. 
The dynamics started from the uniform wavefunction. The moving random potential nucleates 
sound waves as well as vortices, but the incompressible kinetic energy E l kin (t) due to vortices 
becomes dominant in the total kinetic energy Ei^ n (t). The energies soon become almost constant, 
and steady QT is obtained. The energy cascade can be confirmed quantitatively by checking that 
the energy dissipation rate e of E l kin (t) is comparable to the energy flux II(A;, t) via a Richardson 
cascade in the inertial range. The energy flux is found to be approximately independent of k and 
comparable to e. As shown in Fig. 0] (b), the energy spectrum is quantitatively consistent with 
the Kolmogorov law in the inertial range 2tt/Xq < k < 2tt/£, equivalent to 0.79 < k < 6.3. 




The above arguments are chiefly limited to length scales larger than the mean spacing i of a 
vortex tangle, in which the Richardson cascade is effective. What happens at smaller scales where 
the Richardson cascade is less effective? The most probable scenario is the Kelvin-wave cascade. 
A Kelvin wave is a deformation of a vortex line into a helix, with the deformation propagating as 
a wave along the vortex [47j . Kelvin waves were first observed by making torsional oscillations in 
uniformly rotating superfluid 4 He |4g]. The approximate dispersion relation for a rectilinear 
vortex is u>k = (nk 2 )/(4:7r)(ln(l/kao) + c) with a constant c ~ 1 and a cutoff parameter ao equal 
to the vortex core size. Here k is the wavenumber of an excited Kelvin wave, which is different 
than the wavenumber used for the energy spectrum above. At a finite temperature, a significant 
fraction of the normal fluid contains damped Kelvin waves because of mutual friction. At very 
low temperatures, however, mutual friction does not occur and the only possible mechanism of 
dissipation is the radiation of phonons [H(3]. Phonon radiation becomes effective only when the 
frequency is very high, typically of the order of GHz (k ~ 1CP 1 nm _1 ). Therefore, a mechanism 
is required to transfer energy to such high wavenumbers for Kelvin-wave damping. An early 
numerical simulation based on the vortex filament model showed that Kelvin waves are unstable 
to the buildup of sidebands 5l[. Consequently nonlinear interactions between different Kelvin 
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wavenumbers can transfer energy from small to large wavenumbers, which constitutes a Kelvin- 
wave cascade. This idea was first suggested by Svistunov [52(| and later developed and confirmed 
through theoretical and numerical work by Kivotides et al. [53] , Vinen et al. [541 ] , and Kozik and 
Svistunov 55, 5^, 57]. 



Consider the nature of the transition between the Richardson and the Kelvin-wave cascades. 
Turbulence at scales larger than the mean vortex spacing obeys the Kolmogorov law with a 
Richardson cascade. In this regime, the nature of individual vortices is not relevant, and the 
large-scale velocity field created collectively by the vortices mimics classical turbulent behavior. 
This regime may therefore be called classical. On the other hand, at length scales smaller than 
the mean vortex spacing the energy is transferred along each vortex through a Kelvin-wave 
cascade; this regime has no classical analog and can be called quantum. Several theoretical 



considerations of the classical-quantum crossover have been proposed [58|, [59|, [60[] , but this topic 
is not yet settled and remains controversial. 



4-3.2 Quantum turbulence created by vibrating structures 

Recently, vibra ting structures including discs, spheres, grids, and wires, have been used for 
research into QT [61)]. Despite differences between the structures, the experiments show surpris- 
ingly similar behavior. When a classical viscous fluid steadily flows past a structure, the flow 
changes from laminar to turbulence in the wake as the Reynolds number increases. The net drag 
force on the structure is often expressed in terms of the drag coefficient Co as Fr> = 1/2CdpAv 2 , 
where A is the projected area of the structure onto a plane normal to the flow. For laminar vis- 
cous flow, the drag is approximately proportional to v so that Cr> ~ v~ 1 , while turbulent flow 
with high Reynolds number yields a drag coefficient of order unity. This drastic reduction 
in the drag force is called the drag crisis, and is closely related to the formation of the boundary 
layer [62}]. A boundary layer in a high velocity gradient produces eddies. However, such bound- 
ary layers do not appear in the superfluid component. Another classical aspect is the effect of 
oscillations. The oscillatory case is more complicated than that of steady flow because a second 
length scale appears in addition to the linear size of the structure, namely the viscous pene- 
tration depth 5 = \Jlv /u> where to is the oscillation frequency of the structure [6l|. Owing to 
the importance of the topic in marine engineering, a number of experimental studies have been 
published on oscillatory classical flow past a cylinder. As laminar flow becomes unstable, the 
drag coefficient Cd reduces to a value of order unity. However, the values of Cd oscillate within 
the range 0.5 to 2.0, showing a more complicated transition to turbulence than in the case of 
steady flow. Visual observations show that vortices are formed in the wake of the oscillating 
cylinder. A few experimental studies have reported classical flow around an oscillating sphere. A 
typical experiment shows that a single vortex ring is generated and shed during each half period 
of oscillation [i3] • 

Typical behavior appeared in the pioneering observation of QT for an oscillating microsphere 
in superfluid 4 He by Jager et al. [641 ] . The sphere they used had a radius of approximately 
100 /an. It was made from a strongly ferromagnetic material and was magnetically levitated in 
superfluid 4 He. Its response to an alternating drive was observed. At low drives, the velocity 
response v was proportional to the drive Fr>, taking the "laminar" form Fjj = X(T)v with a 
temperature-dependent coefficient A(T). At high drives, the response changed to the "turbulent" 
form Fd = j(T)(v 2 — Vq) above a critical velocity vq. At low temperatures, the transition from 
laminar to turbulent response was accompanied by hysteresis. Subsequently, several groups have 
experimentally investigated the transition to turbulence in superfluid 4 He and 3 He-B by using 
grids, wires, and tuning forks [61]. 

These experimental studies reported some common behaviors independent of the details of 
the structures, including its type, shape, and surface roughness. The observed critical velocities 
were in the range from 1 to 200 mm/s. Since the velocity was lower than the Landau critical 
velocity of approximately 50 m/s, the transition to turbulence should come not from intrinsic 
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nucleation of vortices but from extension or amplification of remnant vortices^ Such behavior is 
seen in the numerical simulation using the vortex filament model [66]. Figure 5 shows how the 
remnant vortices initially attached to a sphere develop into turbulence under an oscillating flow. 
Many unresolved problems remain, such as the nature of the critical velocity and the origin of 
the hysteresis in the transition between the laminar and turbulent response. 





Figure 5. Evolution of a vortex line near a sphere of radius 100 /im in an oscillating superflow of 150 mm/s at 200 Hz. 
[Hanninen, Tsubota, and Vinen, Phys. Rev. B 75 (2007) 064502, reproduced with permission. Copyright 2007 by the 
American Physical Society.] 

4-3.3 Visualization of quantum turbulence 

There has been little direct experimental information about flow in superfluid 4 He because 
standard flow visualization techniques are not applicable to it. However, the situation is changing 
[67| . For QT, one can seed the fluid with tracer particles in order to visualize the flow field. 
Quantized vortices are observable by appropriate optical techniques. 

A significant contribution was made by Zhang and Van Sciver [6^]. Using the particle image 
velocimetry (PIV) technique with 1.7-/im-diameter polymer particles, they visualized a large- 
scale turbulent flow both in front of and behind a cylinder in counterflowing superfluid 4 He at 



1 Superfluid 4 He at rest usually contains remnant vortices [651 . They are metastably pinned between the walls of the container 
and are established when the liquid is cooled down through the A point to the superfluid phase. 
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finite temperatures. In classical fluids, such turbulent structures are seen downstream of objects 
such as cylinders, and the structures periodically detach to form a vortex sheet. In the case of 
4 He counterflow, however, the locations of the large-scale turbulent structures were relatively 
stable; they did not detach and move downstream, although local fluctuations in the turbulence 
were evident. 

Do the tracer particles follow the normal flow, the superflow, or a more complex flow? Barenghi 
et al. studied this question theoretically and numeric ally to show that the situation changes 
depending on the size and mass of the tracer particles foil . 7ol . 71]. 

Another important contribution was the visualization of quantized vortices by Bewley et al. 
[t3 ] . In their experiments, liquid helium was seeded with solid hydrogen particles smaller than 
2.7 /xm in diameter at a temperature slightly above T\, following which the fluid was cooled to 
below T\. When the temperature is above T\, a homogeneous cloud of the particles was seen that 
disperses throughout the fluid. However, on passing through T\, the particles coalesced into web- 
like structures. Bewley et al. suggested that these structures represent decorated quantized vortex 
lines. They reported that the vortex lines form connections rather than remaining separated, 
and were homogeneously distributed throughout the fluid. The observed fork- like structures 
may indicate that several vortices are attached to the same particle. Applying this technique to 
thermal counter flow, Paoletti et al. have directly visualized the two-fluid behavior (73l |. 



5 Quantum turbulence in atomic BECs 

The long history of research into superfluid helium has uncovered two main cooperative phe- 
nomena of quantized vortices: vortex lattices under rotation and vortex tangles in QT. However, 
almost all studies of quantized vortices in atomic BECs have been limited to vortex lattices. 
This section briefly discusses QT in atomic BECs. 

Kobayashi and Tsubota proposed an easy, powerful method to make a steady QT in a trapped 
BEC, by using precession [7J]. The dynamics of the wavefunction are described by the GP 
equation with dissipation. First, one traps a BEC cloud in a weakly elliptical harmonic potential, 

U(r) = ^[(1 - *i)(l - S 2 )x 2 + (1 + - 5 2 )y 2 + (1 + 5 2 )z\ (11) 

where the parameters 5\ and 5 2 exhibit elliptical deformation in the xy- and zx-planes. Sec- 
ond, to transform the BEC into a turbulent state, a rotation is applied along the z-axis fol- 
lowed by a rotation along the shown in Fig. 6(a). The rotation vector is = 
(£l x , sm &xt, cos £l x t) where £l z and Sl x are the frequencies of the first and second rota- 
tions, respectively. Consider the case where the spinning and precessing rotational axes are 
perpendicular to each other. Then the two rotations do not commute and cannot be represented 
by their sum. This form of precession is also used for turbulence in water [751 ] . 

Starting from a stationary solution without rotation or elliptical deformation, Kobayashi et 
al. numerically calculated the time development of the GP equation by turning on the rotation 
£l x = £l z = 0.6 and elliptical deformation 6% = 5 2 = 0.025. By monitoring the kinetic energy 
and anisotropic parameters, the system eventually becomes statistically steady and isotropic. In 
the steady state, the spectrum E r kin (k,t) of the incompressible kinetic energy per unit mass is 
consistent with the Kolmogorov law [Fig. 6(b)]. The inertial range which sustains the Kolmogorov 
law is determined by the Thomas- Fermi radius -Rtf and the coherence length £. The application 
of a combined precession around three axes enables one to obtain more isotropic QT [761 ] . The 
velocity field in a BEC cloud can be observed by Bragg spectroscopy [73], enabling one to obtain 
the energy spectrum. 

There are several advantages of studying QT in atomic BECs compared to in superfluid 
helium. First, one can observe the vortex configuration, probably even the Richardson cascade 
process of vortices. In fact one can directly study the relation between the real-space Richardson 
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Figure 6. QT in atomic BECs. (a) The BEC cloud is subject to precessions, (b) Energy spectrum of a steady QT. The dots 
are a numerically obtained spectrum of the incompressible kinetic energy, while the solid line is the Kolmogorov spectrum. 
Here, Rtf is the Thomas-Fermi radius amd at ra p = \J h/mui is the characteristic length scale of the trap. 



cascade and the wavenumber-space cascade (the Kolmogorov spectrum). Second, one can control 
the transition to turbulence by changing the rotational frequencies or other parameters. For 
example, rotation along one axis forms a vortex lattice. If one applies another rotation, it may 
rotate the lattice if the frequency is low. If the frequency is high, however, the second rotation 
can melt the lattice into a vortex tangle. That would enable one to investigate in detail the entire 
transition to turbulence. Third, by changing the shape of the trapping potential, one can study 
the effect of the anisotropy on turbulence; a typical question is how the Kolmogorov spectrum 
is changed when the BEC becomes anisotropic. This issue is related to the transition between 
2D and 3D turbulence. Fourth, one can create rich QT in multicomponent BECs by controlling 
the interaction parameters. Multicomponent BECs have been extensively studied and found to 
create many kinds of vortex structures depending on the intercomponent interaction (7^]. By 
using such systems, we will be able to make multicomponent QT coupled through the interaction. 



6 Conclusions 

In this article, we have reviewed recent research on quantum hydrodynamics and QT in superfluid 
helium and atomic BECs. Research on QT is currently a key branch in low-temperature physics, 
attracting the attention of many scientists. It is comprised of quantized vortices as definite 
elements, in contrast to conventional turbulence. Thus, the investigation of QT may lead to a 
breakthrough in understanding a mystery of nature pondered since the era of da Vinci. There 
are many related topics not addressed in this article, regarding which we refer the readers to 
other review articles 15. 
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